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Abstract 

In this paper we study the magnetic susceptibihty and other thermodynamic properties of the po- 
larized nuclear matter at finite temperature using the lowest order constrained variational (LOCV) 
method employing the AVi^ potential. Our results show a monotonic behavior for the magnetic 
susceptibility which indicates that the spontaneous transition to the ferromagnetic phase does not 
occur for this system. 
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I. INTRODUCTION 



The magnetic susceptibility is one of the most important magnetic properties of the dense 
matter, its behavior specifies whether the spontaneous phase transition to a ferromagnetic 
state occurs. This transition in nuclear matter could have important consequences for the 
physical origin of the magnetic field of the pulsars which are believed to be rapidly ro- 
;ating neutron stars with strong surface magnetic fields in the range of 10^^ — 10^'^ Gauss 
, Considering different stages of the neutron star formation at different tem- 
peratures, the study of the magnetic properties of the polarized nuclear matter at finite 
temperature is of special interest in the description of protoneutron stars. A protoneutron 
star (newborn neutron star) is born within a short time after the supernovae collapse. In 
this stage the interior temperature of the neutron star matter is of the order 20-50 MeV j^. 
The magnetic susceptibility of nucleonic matter is also an useful quantity to estimate the 
mean free path of the neutrino in the dense nucleonic matter which is relevant information 
for the understanding of the mechanism underlying the supernova explosion and the cooling 



process of the neutron stars 



6|- 



There exists several possibilities of the generation of the magnetic field in a neutron 
star. From the nuclear physics point of view, such a possibility has been studied by several 
authors using different theoretical approaches [7-30], but the results are still contradictory. 
In most calculations, the neutron star matter is approximated by pure neutron matter at 
zero temperature. The properties of the polarized neutron matter both at finite and zero 
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Some calculations show 



30| . Some 



temperature have been studied by several authors [30|, ^ 

that the neutron matter becomes ferromagnetic for some densities {3, S 14, 1^, 
others, using the modern two-body and three-body realistic interactions, show no indication 
of the ferromagnetic transition at any density for the neutron matter and the asymmetrical 



nuclear matter 
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32l |. The results of another calculation show both behaviors; the 



DIP force exhibits a ferromagnetic transition whereas no sign of such transition is found 
for Dl at any density and temperature SjJ- The influence of the finite temperature on the 
antiferromagnetic (AFM) spin ordering in the symmetric nuclear matter with the effective 
Gogny interaction within the framework of a Fermi liquid formalism has been studied by 



Isayev 
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3J]. We note that in the symmetric nuclear matter corresponding to AFM spin 
ordering, we have Ap|| = (p„| + ppi) - {pni + Ppt) 7^ ^"^^ ^PTT = PT ~ Pi = 0- this 
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article, we use the lowest order constrained variational (LOCV) formalism to investigate 
the possibility of the transition to a ferromagnetic phase for the polarized hot symmetrical 
nuclear matter. 



The 



35|,l36|, 



jOCV method has been developed to study the bulk properties of the quantal fluids 



33] . This technique has been used for studying the ground state properties of finite 



nuclei and treatment of isobars 
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40|. Modarres has extended the LOCV method to the 



finite temperature calculations and has applied it to the neutron matter, nuclear matter and 



asymmetrical nuc 



of these systems [41 



ear matter in order to calculate the different thermodynamic properties 



42 



43 



44| . Few years ago, we calculated the properties of nuclear 



matter at zero and finite temperature using the LOCV method with the new nucleon-nucleon 



potentials 



45|, 



46, 



43]. The LOCV method has several advantages with respect to the other 



many-body formalism. These are as follows: (i) Since the method is fully self-consistent, it 
does not introduce any free parameters into the calculations, (ii) It considers the constraint 



in the form of a normalization constraint 
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possible [37 



39 



40 



44 



to keep the higher-order terms as small as 
451] and it also assumes a particular form for the long-range behavior 



of the correlation function in order to perform an exact functional minimization of the two- 
body energy with respect to the short-range behavior of the correlation function, (iii) The 
functional minimization procedure represents an enormous computational simplification over 
the unconstrained methods (i.e. to parameterize the short-range behavior of the correlation 
functions) which attempt to go beyond the lowest order. 



Recently, we have computed the properties of the polarized neutron matter 



49], polar- 



ized symmetrical 50| and asymmetrical nuclear matters 5ll] and also polarized neutron star 



matter 



5l| at zero temperature using the microscopic calculations employing the LOCV 



method with the realistic nucleon-nucleon potentials. We have concluded that the sponta- 
neous phase transition to a ferromagnetic state in these matters does not occur. We have 



also calcu 
perature [ 



ated the thermodynamic properties of the polarized neutron matter at finite tem- 



52] such as the total energy, magnetic susceptibility, entropy and pressure using 



the LOCV method employing the AVis potential [53|. Our calculations do not show any 
transition to a ferromagnetic phase for a hot neutron matter. 

In the present work, we intend to apply the LOCV calculation for the polarized symmet- 
rical nuclear matter at finite temperature using the AVis potential. 



3 



II. FINITE TEMPERATURE CALCULATIONS FOR POLARIZED NUCLEAR 
MATTER WITH THE LOCV METHOD 



We consider a system of A interacting nucleons with A^'^^ spin-up and A^~^ spin-down 
nucleons. For this system, the total number density (p) and spin asymmetry parameter [6) 
are defined as 

P = P^^^+P^-\ 
p(+) _ p{-) 

6 = ^ (1) 

P 

6 shows the spin ordering of the matter which can have a value in the range of 5 = 0.0 (un- 
polarized matter) to 5 = 1.0 (fully polarized matter). To obtain the macroscopic properties 
of this system, we should calculate the total free energy per nucleon, F, 

F = E -TS'^+^ - rS^'\ (2) 

E is total energy per nucleon and 5**^*^ is the entropy per nucleon corresponding to spin 
projection i, 



5«(p,T) = -l5:{[l-n«(fc,r,p«)]ln[l-n»(fc,r,p«)] 

^ k 

+n^'\k, T, p«)lnn«(A;, T, p«)}. (3) 



where n^^\k,T, p^^^) is the Fermi-Dirac distribution function, 

n^'\k,T,p^'^) = ^^[,(i)(fc,r,pW)-M«(r,pW)] i' (^) 

In the above equation l3 = , p^*^ being the chemical potential which is determined at 
any adopted value of the temperature T, number density p'-*-' and spin polarization S, by 
applying the following constraint, 

5:n«(fc,r,p«)=A«, (5) 

k 

and e*^*-* is the single particle energy of a nucleon. In our formalism, the single particle energy 
of a nucleon with momentum k and spin projection i is approximately written in terms of 



the effective mass as follows 
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34| 



^^^Hfe,T,p(-)) = ^J^^^^^^ +^(-)(T,p(-)). (6) 



In fact, we use a quadratic approximation for single particle potential incorporated in the 
single particle energy as a momentum independent effective mass. f/^*^(T, p*^*)) is the mo- 
mentum independent single particle potential. We introduce the effective masses, Tn*^^\ as 



variational parameters 
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541 ]. We minimize the free energy with respect to the varia- 
tions in the effective masses and then we obtain the chemical potentials and the effective 
masses of the spin-up and spin-down nucleons at the minimum point of the free energy. This 
minimization is done numerically. 

As it is also mentioned in the pervious section, for calculating the total energy of the 
polarized symmetrical nuclear matter, we use the LOCV method. We adopt a trial many- 
body wave function of the form 

^ = (7) 

where (j) is the uncorrelated ground state wave function (simply the Slater determinant of 
plane waves) of A independent nucleons and = J^[l---A) is an appropriate A-body 
correlation operator which can be replaced by a Jastrow form i.e., 

J^ = Sl[f{zj), (8) 

i>j 

in which 5 is a symmetrizing operator. Now, we consider the cluster expansion of the energy 
functional up to the two-body term jol. 

For the hot nuclear matter, the one-body term Ei is 

^1 = ^;+)+^;-), (10) 



where 



The two-body energy E2 is 



where 



E? = Y.^n(^{k,T,p(^^). (11) 



^2 = ^Y.i'jHl2)\tj-jt), (12) 



K12) = -|^[/(12), [V?2, /(12)]] + /(12)V^(12)/(12). (13) 



In above equation, /(12) and V{12) are the two-body correlation and potential. In our 



calculations, we use the AVis two-body potential which has the following form 



53|, 



18 



V{l2) = J2v^^\n2)0 
p=l 



(p) 

12 5 



(14) 



where 



^(p=l-18) 



1, (Ti • (72, Ti • r2, ((Ti • (T2) (ti • T-2) , S12, S'i2(ri • T-^) , 

L ■ S, L ■ S(ri ■ r2), L^ L2(ai • ^2), L2(ri ■ r2), 
L^K-^2)(ri-r2), (L■S)^ (L ■ S)2(ri ■ r2), 

T12, (O"! ■ 0-2)Ti2, 512X12, {t^1+T^2)- 



(15) 



In above equation, 5*12 = [3(cri-f)((T2-'r)— (Ti-(T2] is the tensor operator and T12 = [3(ri-f)(r2- 
r)—Ti-T'2\ is the isotensor operator. The above 18 components of the AVi^ two-body potential 
are denoted by the labels c, a, r, ar, t, tr, /s, /sr, /2, /2cr, /2r, /2crr, /s2, Zs2r, T, crT, tT and r^r, 
respectively [ssj. In the LOCV formalism, the two-body correlation /(12) is considered as 
the following form 37l |. 



/(12) 



E/^'Hn2)Pi?, 
fc=i 



(16) 



where 



p(fc=l-3) 

-'12 



1 1 

4 ~ 4 
1 1 

4 ^ 12 



-Oyi ] ) 



12 ) 



r)(5) 
•-^12 



(17) 



The operators Oy^ and Oy^ are given in Eq. ( fTSl) . Using the above two-body correlation 
and potential, after doing some algebra we find the following equation for the two-body 
energy, 

J_ (fl\ ^ (2J + 1)(2T+1^ 
^'pWI jLk.T. 2(25 + 1) 



■[l-(-l) 



L+S+Ti 



11 ^2 

22' 



X dr 



f^^"a^l^\kjr) + 1^ ({K - 3K + (K - 3K.)(4T - 3) 
HVt - 3Kr)(4T)}a«'( V) + [V12 - 



Wnr - 3l^2..)(4T - 3)]c« ik^r)) {f^^^f 

fc=2,3 



6A; + 14)V, + -(A;-1)V^, 



+[K + V,r + (-6A; + U)Vu -{k- l)l^.r](4T - 3) 
+[Vt + V^r + (-6fc + U)VtT] [T{QT^ - 4)] + 

+ [Vi2 + Vi2a + iVi^r + V^2.r)(4T - 3)]c«'(V) 

+[(1^.2 + T^.20(4T-3)]d«'(V))/«'; 

+ -^[[(V,.. - 2iyi2ar + - 3Vz,2r)(4T - 3)] 

where a = {J, L, 5", -S^} and the coefficient a^^^^, etc., are as follows, 

af\x) - x2[7/j_i,5.(x) 

- X^[/323/j-l,5.(x) - /323/j+l,5.(x)], 

cW^(x) = x^vJl^sA^), 
ci')'(x) = a;2[772/j_i,s.(x) + V2lj+i,sA^)l 

c(f)'(x) = X2[7y3/j-l,5.(2:) + 1^3/7+1,5.(2^)], 
df\x) = X^[6/J-1,5.(X) + A2/J+1,5.(X)], 
d(f)'(x) = X2[e3/j-l,5.(x) + A3/J+1,5.(X)]. 

In above equations, we have 

2J+1' ^ 2J+1' 2J+1 ' 



(29) 



V2 



J{J^ + 2J+1) 
2J + 1 



J( + J + 2) 
2J+ 1 



(30) 



J3 + 2 J2 + 2 J + 1 J( J2 + J + 4) 

« — r, T , 1 ' S3 — 



2J+ 1 



2J + 1 



(31) 



_J(J2 + J+1) _ j3_^2j2 + 5J + 4 



2J+ 1 



2J + 1 



(32) 



and 



1 



k\dkikldk2ni{ki,T , pi)nj{k2,T, pj)Jj{\k2 - ki\r), 



(33) 



27rV 

where Jj{x) is the Bessel's function . 

Now, we minimize the two-body energy Eq. f|T8|) with respect to the variations in the 
correlation functions fQ^''\ but subject to the normahzation constraint jsT, 46|, 

1 



A 



E(^J>1-/'(12) ^J)a = 0■ 



(34) 



^3 



In the case of polarized symmetrical nuclear matter, the Pauli function hs^{r) is as follows 



hsjr) 



'1/2 



; S, = ±l 



(35) 



where 



7»(r) 



(36) 



From the minimization of the two-body cluster energy we get a set of coupled and uncoupled 
Euler-Lagrange differential equations. The Euler-Lagrange equations for uncoupled states 



are 



(1)" 

f ag 



m 
1? 



■[K-3K + (K-3V;.)(4T-3) 
^(Vt - 3\4t)[T(6T2 - 4)] + 2K.T, + A] 

+ f (Vi2 - 3^,2. + {Vl2r - 3l^2.r)(4T - 3))4^}# = 0, 



(37) 
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while for the coupled states, these equations are written as follows, 



9, 



a 



{^ + §[V, + V, + 2Vt- Vis + (K + V^r + 2Vtr - 
+ {Vt + + 2VtT)[T{6T^ - 4)] + + A] 

+ §[Vl2 + Vl2a + {Vl2r + V/2<xr) (4T - 3)]^ 
+ §[Vls2 + V1S2MT -3)]'^ + 



l^,,)(4T-3) 



+ {yisr - 2Vl2r - 2Vnar ' 



3V^.2.)(4T-3)]}^<7, 



(3) 







(38) 



- + f [V; + K - 4Vi - 2Vu + (K + - 4Ft, - 2Vur)(^T - 3) 

+(Vt + - 4ViT) [r(6T2 - 4)] + 2v;,r, + A] 

+ f [^2 + Vl2. + (V^/2r + V52..)(4T - 3)]^ 

+f [^s2 + v^.2.(4T - 3)]4;i + 

+{^ - #[V^s - IVn - 2Vi2. - Wis2 

HVlsr - 1Vl2r - 1Vl2ar ' 3^.2.) (4T - 3)]}^#^# = 0, (39) 



The primes in the above equations mean differentiation with respect to r. The Lagrange 
multiplier A is introduced by the normalization constraint, Eq. flM|) . Now, we can calculate 
the correlation functions by numerically solving these differential equations and then using 
these correlation functions, the two-body energy is obtained. Finally, we can compute the 
energy and then the free energy of the system. 

III. RESULTS AND DISCUSSION 

We have presented the effective masses of the spin-up and spin-down nucleons as functions 
of the spin polarization (5) at p = O.bfm"^ and T = 20 MeV in Fig. [H It is seen that 
the difference between the effective masses of spin-up and spin-down nucleons increases by 
increasing the polarization. We also see that the effective mass of spin-up nucleons increases 



where 



9, 



(40) 
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by increasing the polarization whereas the effective mass of spin-down nucleons decreases 
by increasing the polarization. These behaviors have been also seen for the effective mass 



of the neutron in the case of spin polarized hot neutron matter 3(J, |3l|, l32|, |52| . A similar 



qualitative behavior of the nucleon effective mass as a function of the isospin asymmetric 
parameter j3 = has been already found in non-polarized isospin asymmetric nuclear 



matter, as it has been already discussed in Refs. [55|, |56 |. 

The free energy per nucleon of the polarized hot nuclear matter versus the total number 
density (p) for different values of the spin polarization (6) at T = 10 and 20 MeV is shown 
in Fig. [21 It can be seen that for each value of the temperature, the free energy increases by 
increasing both density and polarization. From Fig. [2], it is seen that the free energy of the 
polarized hot nuclear matter decreases by increasing the temperature. We have seen that 
above a certain values of the temperature and spin polarization, the free energy does not 
show any bound states for the polarized hot nuclear matter. From Fig. [2] we also see that for 
all given temperatures there is no crossing of the free energy curves for different polarizations 
and the difference between the free energy of the nuclear matter at different polarizations 
increases by increasing the density. This indicates that the spontaneous transition to the 
ferromagnetic phase does not occur in the hot nuclear matter. We have also compared the 
free energy per nucleon of the unpolarized case of the nuclear matter at different tempera- 
tures in Fig. [3l It is seen that the free energy of unpolarized nuclear matter decreases by 
increasing the temperature. We can see that above a certain temperature, the free energy 
does not show the bound state (the minimum point of the free energy) for the unpolarized 
nuclear matter. 

For the polarized hot nuclear matter, the magnetic susceptibility, x, which characterizes 
the response of the system to the magnetic field can be calculated by using the following 
relation, 

X = 7^. (41) 



.dP J 5=0 

where /i is the magnetic moment of the nucleons. Fig. S] shows the ratio xf/x ^is a function 
of the temperature at p = 0.16/m~'^, where xf is the magnetic susceptibility for a nonin- 
teracting Fermi gas. As it can be seen from this figure, this ratio is inversely proportional 
to absolute temperature without any anomalous change in its behavior. This indicates that 
hot nuclear matter is paramagnetic. The ratio Xf/x has been also shown versus the to- 

10 



tal number density at temperature T = 20 MeV in Fig. [51 A magnetic instability would 
require X-f/x < 0- It is seen that the value of Xf/x is always positive and monotonically 
increasing up to highest density and does not show any spontaneous phase transition to the 
ferromagnetic phase for the hot nuclear matter. 

The difference between the entropy per nucleon of the fully polarized and unpolarized 
cases of the nuclear matter is plotted as a function of the total number density at T = 20 
MeV in Fig. [61 It is seen that for all given values of the density, this difference is negative. 
This shows that the fully polarized case of hot nuclear matter is more ordered than the 
unpolarized case. We also see that the magnitude of this deference decreases by increasing 
the density. The entropy per nucleon of the polarized hot nuclear matter versus the spin 
polarization for fixed density p = 0.5/m~^ and temperature T = 20 has been presented 
in Fig. [3 It is shown that the entropy decreases by increasing the polarization. It is also 
shown that the highest value of the entropy occurs for the unpolarized case of the hot nuclear 
matter. For the polarized hot nuclear matter, the following condition for the effective mass 
prevents the anomalous behavior of the entropy versus the spin polarization 30j, 

m*{p,5= 1.0) 



m*{p,6 = 0.0) 



< 2^/3, (42) 



where m*{p,S = 1.0) and m*{p,6 = 0.0) are the effective masses of the fully polarized 
and unpolarized nuclear matter, respectively. This condition was first derived in Ref. [s^ 
for the particular case of the Skyrme interaction where the effective mass is independent 
of the momentum and temperature and therefore the single particle potential is purely 
parabolic. In our approach, the effective mass depends on both density and temperature 
but is independent of the momentum. In other words, a similar rigorous condition can not be 
obtained straightforwardly. However, within this approximation, one can use the condition 
of Eq. (H2l) . From our result for the effective mass at T = 20MeV for p = O.S/m"^ (Fig- 
[1]), we have found that this ratio is 1.24. We see that this value is smaller than the above 
limiting value which indicates that the entropy of polarized case of the hot nuclear matter is 
always smaller than the entropy of unpolarized case. This so-considered "natural" behavior 
was also found in the case of Gogny 3l|| and in the BHF analysis of Ref. 32|. In contrast, 
for Skyrme forces the entropy per particle of the polarized phase is seen to be higher than 
the non-polarized one above a certain density jsi]. 

Finally, we have plotted the pressure of the polarized hot nuclear matter as a function 
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of the total number density (p) for different polarizations at T = 10 and 20 MeV in Fig. 
[HI For all values of temperature and polarization, it is seen that the pressure increases by 
increasing the density. For this system, we see that at each temperature the equation of 
state becomes stiffer as the polarization increases. For each polarization, it is found that 
the pressure of the polarized hot nuclear matter increases by increasing the temperature. 

IV. SUMMARY AND CONCLUSIONS 

The lowest order constrained variational (LOCV) method has been used for calculating 
the susceptibility of the polarized hot nuclear matter and some of the thermodynamic prop- 
erties of this system such as the effective mass, free energy, entropy and the equation of 
state. In our calculations, we have employed the AVis potential. Our results show that the 
spontaneous transition to the ferromagnetic phase does not occur for the hot nuclear matter. 
We have seen that the spin polarization substantially affects the thermodynamic properties 
of the hot nuclear matter. 
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FIG. 1: The effective mass of spin-up (full curve) and spin-down (dashed curve) nucleons versus 
the spin polarization (S) for density p = 0.5/m~^ at T = 20 MeV. 
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FIG. 2: The free energy per nuclcon of the polarized hot nuclear matter as a function of the total 
number density (p) for different values of the spin polarization {6) at 7" = 10 (a) and T = 20 MeV 
(b). 
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FIG. 3: The free energy per nucleon of the nuclear matter versus the total number density (p) for 
unpolarized case at T = 0, 10 and 20 MeV . 




FIG. 4: The magnetic susceptibility of the hot nuclear matter versus the temperature at p 
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FIG. 5: The magnetic susceptibility of the hot nuclear matter versus the total number density (p) 
at T = 20 MeV. 
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FIG. 6: As Fig. 4 but for the entropy difference of the fully polarized and the unpolarized cases. 
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FIG. 7: The entropy per nucleon as a function of the spin polarization {S) for density p = 0.5/m 
at r = 20 MeV . 
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FIG. 8: The equation of state of the polarized hot nuclear matter for different values of the spin 
polarization (S) at T = 10 (a) and T = 20 MeV (b). 
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